A. 
A. 
N. 



H 



Centrality dependence of charged hadron production 
in deuteron+gold and nucleon+gold collisions at ^/snn = 200 GeV 

S.S. Adler,5 S. Afanasiev,20 C. Aidala,i° N.N. Ajitanand,44 Y. Akiba,2i'40 a. Al-Jamel,^^ J. Alexander,44 
K. Aoki,^'^ L. Aphecetche,''^ R. Armendariz,^^ S.H. Aronson,^ R. Averbeck,''^ T.C. Awes,^^ V. Babintsev,^"^ 
Baldisseri," K.N. Barish,^ P.D. Barnes,^^ B. Bassalleck,^^ S. Bathe,^'^! S. Batsouli,!'^ V. Baublis,^^ F. Bauer,^ 
Bazilevsky,5'4i S. Belikov,!^' "^B M.T. Bjorndal,i" J.G. Boissevain,^^ H. Borel," M.L. Brooks,^^ D.S. Brown,^^ 
Bruner,34 D. Bucher,^! H. Buesching,^'^! V. Bumazhnov," G. Bunce,^''*! J.M. Burward-Hoy,^*'^'^ S. Butsyk,'^^ 
X. Caniard,'^^ Chand,"* W.C. Chang,^ S. Chcrnichcnko,!^ C.Y. Ghi,^° J. Chiba,^! M. Chiu,i° I.J. Choi,53 

R.K. Choudhury 4 T. Chujo,^ V. Cianciolo,^^ Z. Citron ,^5 Y. Cobigo," B.A. Cole/o M.P. Comets,^^ 
Constantin,!^ M. Csanad,!^ t. Csorgo,^^ j Cussonneau,'^^ d'Enterria,!" K. Das," G. David,^ F. Deak," 
Delagrange,^^ A. Denisov,^^ A. Deshpande,^^ E.J. Desmond,^ A. Devismes,''^ O. Dietzsch,''^ J.L. Drachenberg/ 

O. Drapier,26 A. Drees,'^^ A. Durum,i^ D. Dutta,"* V. Dzhordzhadze,"*^ Y.V. Efremenko,^^ H. En'yo,'*"''*^ 
B. Espagnon,37 S. Esumi^^ D.E. Fields,34^4i q Finck,*^ F. Fleuret,^^ S.L. Fokin,^^ B.D. Fox^i Z. Fraenkel,^^ 
J.E. Frantz,io A. Franz,^ A.D. Frawley,!'' Y. Fukao,25'40,4i g _y. Fung,^ S. Gadrat,^^ M. Germain^e 
A. Glenn 47 M. Gonin,^^ J. Gosset," Y. Goto,'*"'''! R. Granier de Cassagnac,^^ N. Grau,!^ S.V. Greene,^^ 
M. Grosse Perdekamp/^-^'i H.-A. Gustafsson,30 T. Hachiya,!^ J.S. Haggerty,^ H. Hamagaki,* A.G. Hansen,^ 
.P. Hartouni,27 M. Harvey,^ K. Hasuko,'*" R. Hayano,^ X. Re,^^ M. Heffner,^^ T.K. Hemniick,''^ J.M. Heusei 
P. Hidas,22 H. Hiejima,!^ J.G. Hill,!^ R. Hobbs,^'* W. Holzmann,44 K. Homma,!^ B. Hong,^^ A PTnn.mr 35 
Horaguchi 40.41:48 t. Ichihara 40'4i y.V. Ikonnikov,24 K. Imai,25.40 m. Inaba49 M. Inuzuka,* 
Isenhower/ M. Ishihara40 M. Issah,44 A. Isupov,20 B.V. Jacak 45.0 J. Jia45 Q. Jinnouchi,4".4i 

r. ^ T 1 97 n T ..^9 t :\7 ■■^ . S r. . . . • S ."SI . -1 . . 4-0 4-S n 



, 50 

M. Grosse Perdekamp,i8,4i h.-A. Gustafsson,30 T. Hachiya,!^ J.S. Haggerty,^ H. Hamagaki,^ A.G. Hansen,^^ 
E.P. Hartouni,27 M. Harvey,^ K. Hasuko,40 R. Hayano,^ X. He,!^ M. Heffner,^^ T.K. Hemniick,45 J.M. Heuser 4° 

P. Hidas,22 H. Hiejima,!^ J.G. Hill,!^ R. Hobbs,34 W. Holzmann,44 K. Homma,!^ B. Hong,^^ A. Hoover,^^ 
T. Horaguchi 40.41.48 t. Ichihara 40.4i y.V. Ikonnikov,24 K. Imai,25.40 m. Inaba49 M. Inuzuka,* D. Isenhower,i 
L. Isenhower/ M. Ishihara40 M. Issah,44 A. Isupov,20 B.V. Jacak 45.0 J. Jia45 Q. Jinnouchi,40.4i b.M. Johnson,5 
S.G. Johnson,27 K.S. Joo,^^ D. Jouan,^^ F. Kajihara,^ S. Kametani,^. 5i N. Kamihara 40. 48 m. Kaneta,4i 
J.H. Kang,53 K. Katou,5i T. Kawabata,^ A.V. Kazantsev,24 S. Kelly,^.^o ^ Khachaturov,52 A. Khanzadeev,^^ 
J. Kikuchi,5i D.J. Kim,53 E. Kim 43 G.-B. Kim,^^ H.J. Kim,53 E. Kinney,^ A. Kiss,i3 E. Kistenev,5 A. Kiyomichi,40 

C. Klein-Boesing,3i H. Kobayashi,4i L. Kochenda,39 V. Kochetkov,^'' R. Kohara,^^ B. Komkov,3^ M. Konno 4^ 
D. Kotchetkov,^ A. Kozlov,52 p.j. Kroon,5 C.H. Kuberg/'B G.J. Kunde,^^ K. Kurita 40 M.J. Kweon,23 Y. Kwon,53 



M.J 

A 

M.G 



lem-Boesmg, H. Kobayashi, L. Kochenda, V. Kochetkov, R. Kohara, B. Komkov, M. Konno, 
chetkov,^ A. Kozlov,52 P.J. Kroon,5 C.H. Kuberg,i.B G.J. Kunde,^^ K. Kurita 40 M.J. Kweon,23 Y. Kwon,5 
G.S. Kyle,35 R. Lacey 44 J.G. Lajoie,!^ Y. Le Bornec,37 A. Lebedev/^' ^4 S. Leckey 45 D.M. Lee,^^ 
J. Leitch,28 M.A.L. Leite,42 X.H. Li,^ H. Lim 43 A. Litvinenko,20 M.X. Liu,^^ C.F. Maguire,50 Y.I. Makdisi,5 
L Malakhov,20 V.I. Manko,24 Y. Mao,3«.40 G. Martinez,46 H. Masui49 F. Matathias,45 T. Matsumoto,^' 5i 
[.C. McCain,! P.L. McGaughey,^^ Y. Miake,49 T.E. Miller,50 A. Milov,45 S. Mioduszewski,5 G.C. Mishra,i5 
J.T. Mitchell, 5 A.K. Mohanty,4 D.P. Morrison,5 J.M. Moss,^^ D. Mukhopadhyay,52 M. Muniruzzaman,'^ 
S. Nagamiya,2i J.L. Nagle,^. ^o T. Nakamura,!^ J. Newby,47 A.S. Nyanin,24 J. Nystrand,30 E. 0'Brien,5 
C.A. Ogilvie,!^ H. Ohnishi,40 I.D. Ojha,3.50 h. Okada,25.40 k. Okada,40.4i A. Oskarsson,30 I. Otterlund,30 
K. Oyama,^ K. Ozawa,^ D. Pal,52 a.P.T. Palounek,^^ V. Pantuev,45 V. Papavassihou,35 J. Park,43 W.J. Park,23 
S.F. Pate,35 H. Pei,!^ V. Penev,^^ J.-C. Peng,!^ H. Pereira," V. Peresedov,20 A. Pierson,34 C. Pinkenburg,5 
R.P. Pisani,5 M.L. Purschke,5 A.K. Purwar,45 J.M. Quails, ^ J. Rak,!^ I. Ravinovich,52 K.F. Read,36.47 m. Reuter,45 
K. Reygers,3i V. Riabov,39 Y. Riabov,39 G. Roche,^^ A. Romana,26.B M. Rosati,i^ S.S.E. Rosendahl,30 p. Rosnet,^^ 

V.L. Rykov,40 S.S. Ryu,53 N. Saito,25. 40, 4i Sakaguchi,^. 5i S. Sakai,49 V. Samsonov,39 L. Sanfratello,34 
R. Santo,3i H.D. Sato,25.4o s. Sato,5.49 s. Sawada,^! Y. Schutz,46 V. Semenov,!^ R. Seto,^ T.K. Shea,5 1. Shein," 
T.-A. Shibata,40.48 K. Shigaki,!^ M. Shimoniura,49 A. Sickles,45 C.L. Silva,42 D. Silvermyr,^^ K.S. Sim,23 
A. Soldatov,!^ R.A. Soltz,^^ W.E. Sondheim,^^ S.P. Sorensen,47 LV. Sourikova,5 F. Staley," P.W. Stankus,36 
E. Stenlund,30 M. Stepanov,35 A. Ster,^^ S.P. Stoll,5 T. Sugitate,!^ J.P. Sullivan,^^ S. Takagi,49 E.M. Takagui,42 
A. Taketani,40.4i k.H. Tanaka,^! Y. Tanaka,33 K. Tanida,40 M.J. Tannenbaum,5 A. Taranenko,44 p. Tarjan,!^ 
T.L. Thomas,34 Togawa,25.40 j. Tojo,40 H. Torii,25.4i r.s. Towell,^ V-N. Trani,^^ L Tserruya,^^ 
Y. Tsuchimoto,!^ H. Tydesj6,30 N. Tyurin,!'^ T.J. Uam,32 j Vclkovska,5 M. Velkovsky,45 V. Veszprenii,!^ 
A.A. Vinogradov,24 M.A. Volkov,24 E. Vznuzdaev,39 X.R. Wang,i5 Y. Watanabe,40.4i s.N. White,5 
N. Willis,37 F.K. Wohn,i9 C.L. Woody,5 W. Xie,^ A. Yanovich," S. Yokkaichi,40. 4i Q.R. Young,36 
LE. Yushmanov,24 W.A. Zajc,io C. Zhang,io S. Zhou,^ J. Zimanyi,22.B L. Zolin,20 X. Zong,!^ and H.W. vanHecke^^ 

(PHENIX Collaboration) 

^Abilene Christian University, Abilene, TX 79699, USA 
"^Institute of Physics, Academia Sinica, Taipei 11529, Taiwan 
^Department of Physics, Banaras Hindu University, Varanasi 221005, India 
^'Bhabha Atomic Research Centre, Bombay 400 085, India 
Brookhaven National Laboratory, Upton, NY 11973-5000, USA 



2 



^University of California - Riverside, Riverside, CA 92521, USA 
''china Institute of Atomic Energy (CIAE), Beijing, People's Republic of China 
^ Center for Nuclear Study, Graduate School of Science, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 

'^University of Colorado, Boulder, CO 80309, USA 
^"Columbia University, New York, NY 10027 and Nevis Laboratories, Irvington, NY 10533, USA 
" Dapnia, CEA Saclay, F-91191, Cif-sur-Yvette, France 
^'"^ Debrecen University, H-4010 Debrecen, Egyetem ter 1, Hungary 
^^ELTE, Edtvos Lordnd University, H - 1117 Budapest, Pdzmdny P. s. 1/A, Hungary 
Florida State University, Tallahassee, FL 32306, USA 
Georgia State University, Atlanta, CA 30303, USA 
Hiroshima University, Kagamiyama, Higashi- Hiroshima 739-8526, Japan 
^ 'IHEP Protvmo, State Research Center of Russian Federation, Institute for High Energy Physics, Protvino, 142281, Russia 
''^University of Illinois at Urbana- Champaign, Urbana, IL 61801, USA 
'^'^ Iowa State University, Ames, I A 50011, USA 
^" Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 
^'KEK, High Energy Accelerator Research Organization, Tsukuba, Ibaraki 305-0801, Japan 
^^KFKI Research Institute for Particle and Nuclear Physics of the Hungarian Academy 
of Sciences (MTA KFKI RMKI), H-1525 Budapest 114, POBox 49, Budapest, Hungary 
^''Korea University, Seoul, 136-701, Korea 
'^'^ Russian Research Center "Kurchatov Institute", Moscow, Russia 
^^Kyoto University, Kyoto 606-8502, Japan 
Laboratoire Leprince-Ringuet, Ecole Poly technique, CNRS-IN2P3, Route de Saclay, F-91128, Palaiseau, France 
^''Lawrence Livermore National Laboratory, Livermore, CA 94550, USA 
"^^Los Alamos National Laboratory, Los Alamos, NM 87545, USA 
^^LPC, Universite Blaise Pascal, CNRS-IN2P3, Clermont-Fd, 63177 Aubiere Cedex, France 
'^"Department of Physics, Lund University, Box 118, SE-221 00 Lund, Sweden 
Institut fiir Kemphysik, University of Muenster, D-48149 Muenster, Germany 
'^^Myongji University, Yongin, Kyonggido 449-728, Korea 
'^''Nagasaki Institute of Applied Science, Nagasaki-shi, Nagasaki 851-0193, Japan 
University of New Mexico, Albuquerque, NM 87131, USA 
^^New Mexico State University, Las Cruces, NM 88003, USA 
•'^Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA 
'^''IPN-Orsay, Universite Pans Sud, CNRS-IN2P3, BPl, F-91406, Orsay, France 
Peking University, Beijing, People's Republic of China 
'^^PNPI, Petersburg Nuclear Physics Institute, Gatchina, Leningrad region, 188300, Russia 
^°RIKEN (The Institute of Physical and Chemical Research), Wako, Saitama 351-0198, JAPAN 
^'RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973-5000, USA 
'''^ Universidade de S ao Paulo, Instituto de Fisica, Caixa Postal 66318, S ao Paulo CEP05315-970, Brazil 
^'^ System Electronics Laboratory, Seoul National University, Seoul, South Korea 
Chemistry Department, Stony Brook University, Stony Brook, SUNY, NY 11794-3400, USA 
Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, NY 11794, USA 
^'^SUBATECH (Ecole des Mines de Nantes, CNRS-IN2P3, Universite de Nantes) BP 20722 - 44307, Nantes, France 

^''University of Tennessee, Knoxville, TN 37996, USA 
Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo 152-8551, Japan 
Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
''" Vanderbilt University, Nashville, TN 37235, USA 
Waseda University, Advanced Research Institute for Science and 
Engineering, 17 Kikui-cho, Shmjuku-ku, Tokyo 162-0044: Japan 
^^Weizmann Institute, Rehovot 76100, Israel 
^^Yonsei University, IPAP, Seoul 120-749, Korea 
(Dated: June 12, 2008) 

We present transverse momentum {pt) spectra of charged hadrons measured in deuteron-gold 
and nucleon-gold collisions at ^s„„ = 200 GeV for four centrality classes. Nucleon-gold collisions 
were selected by tagging events in which a spectator nucleon was observed in one of two forward 
rapidity detectors. The spectra and yields were investigated as a function of the number of binary 
nucleon-nucleon collisions, i/, suffered by deuteron nucleons. A comparison of charged particle yields 
to those in p+p collisions show that yield per nucleon-nucleon collision saturates with u for high 
momentum particles. We also present the charged hadron to neutral pion ratios as a function of pr- 
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I. INTRODUCTION 

The measurement of hadron spectra at high transverse 
momenta (pt) has been a powerful tool in studying the 
modification of particle production and propagation in 
the nuclear medium at the BNL Relativistic Heavy Ion 
Collider (RHIC). One of the most intriguing results ob- 
served so far at RHIC is the strong suppression of the 
yield of hadrons with px above ~ 2 GeV/c in midcentral 
and central Au+Au collisions relative to the correspond- 
ing yield in p + p collisions scaled by the sum of inde- 
pendent nucleon-nucleon collisions in the Au-|-Au inter- 
action P, '2', . This is consistent with a picture in which 
the hard scattered parton loses significant energy when 
it traverses the medium created during the collision 
There is additional evidence for the strong interaction 
of partons with the final state medium. For example, 
two particle azimuthal correlations demonstrate signifi- 
cant modifications of jets such as width and particle con- 
tent (sj, Q- In the Pt region, 2 < px < 5GeV/c, an 
anomalously large baryon/meson ratio is observed 
The coalescence of thermal and shower partons offers a 
possible explanation of this phenomenon as a final state 
effect 0. 

The hadron spectrum is sensitive to final state effects 
including jet quenching, parton recombination and scat- 
tering of produced particles, in addition to initial state ef- 
fects, such as gluon saturation and nuclear shadowing Q, 
and the Cronin effect [l^, [ll| . If shadowing or gluon sat- 
uration are actually responsible for the large suppression 
seen in Au-|-Au collisions, these same effects would also 
produce a suppression of 20-30% in the yield of high pT [1| 
hadrons even in d-|-Au collisions. 

To study the role of initial state effects in the high 
Pt suppression, d-l-Au collisions at RHIC serve as a con- 
trol experiment, since no dense matter is expected to be 
formed in these collisions. Earlier measurements on in- 
clusive hadron production in d-|-Au collisions show no 
suppression of high pT hadrons [H, [l^. Rather, there 
is a small enhancement in the yield of high momentum 
particles in d-l-Au collisions. This result implies that the 
effects seen in central Au-|-Au collisions arc largely due 
to final state interactions in the dense medium. 

In this paper, we extend our investigation of charged 
hadron production in d+Au collisions at y/s^ — 200 
GeV at midrapidity with detailed centrality selections 
and with significantly higher statistics. The centrality se- 
lection allows one to make quantitative statements about 
the impact parameter dependence of particle production. 

The nuclear modification of high pT hadron produc- 
tion in proton reactions with heavy nuclei has previo usly 
been studied at lower center of mass energies [13, [Hf- 
It was found that hadron yield at high pT in p+A col- 
lisions increases faster than the nuclear mass, A. This 
effect, known as the Cronin effect, is conventionally 
parametrized as A". At y/s ^10-30 GeV, a depends 



on Pt and is found to be greater than unity for pt > 
2GeV/c. The exponent a also depends on the parti- 
cle species; it was found to be larger for baryons than for 
mesons. The Cronin effect is often attributed to the mul- 
tiple scattering of projectile partons propagating through 
the target nucleus |14| . 

A deuteron projectile has an additional interesting 
characteristic: since the average separation between the 
proton and neutron in the deuteron is large, on the order 
of few a fm, there is a class of events in which only one of 
the two nucleons participates in the collision. Thus d-t-Au 
events can be categorized into "p"-|-Au or "n"-|-Au col- 
lisions, when only one of the constituent nucleons inter- 
acts, and pn+Au, when both nucleons participate in the 
collision. In such events the interacting proton or neutron 
is not a free proton or neutron, so it is not exactly the 
same as a p+Au or n+Au interaction. For simplicity we 
will refer to the tagged samples as p-|-Au and n-|-Au and 
collectively 7V-|-Au. The selection of p-|-Au and n-f Au 
collisions in the PHENIX experiment is done by tagging 
events where the incoming deuteron has a non-interacting 
neutron or proton. Combined with centrality selection, 
such event tagging provides better control of the colli- 
sion geometry and of the number of subsequent inelastic 
nucleon-nucleon interactions in the gold nucleus. 

Very little experimental information has been pub- 
lished about neutron-nucleus collisions at high energy. 
Tagging allows the direct comparison of pn+Au, p+Au, 
and n+Au interactions. Results from these collisions are 
compared with p -I- p data measured by the PHENIX ex- 
periment at the same beam energy to shed light on the 
Cronin enhancement and other nuclear effects. Further, 
we have investigated the Cronin effect as a function of 
the number of collisions per participant nucleon in the 
deuteron. To estimate the fraction of pions in the charged 
hadron spectrum, we have looked at hadron to pion ratios 
as a function of transverse momentum. 

The paper is organized as follows. Section II gives a de- 
tailed account of the analysis in which we describe the de- 
tector, centrality selection, charged particle background 
rejection, spectra corrections and systematic errors. In 
Sec. Ill, the centrality and px dependence of the charged 
hadron spectra, and a comparison with the p + p and the 
n'^ data are discussed. A summary is given in Sec. IV. 



II. EXPERIMENT AND ANALYSIS 
A. Detector 

The PHENIX detector consists of four spectrometer 
arms positioned around the vertex of colliding deuteron 
and gold nuclei — two central arms at midrapidity and 
two muon arms at forward rapidities — and a set of de- 
tectors. A detailed description of the detector can be 
found elsewhere, see Ref. [l^ and references therein. 
This analysis uses the two central arms and the set of 
global detectors. 
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Each central arm covers the pseudorapidity range \r]\ < 
0.35, and 90° in azimuth. Charged particles are tracked 
by drift chambers (DCs) located 2 m from the vertex and 
layers of pad chambers (PCI and PCS) positioned at 2.5 
and 5 m in the radial direction. The central spectrometer 
provides axial magnetic fields along the beam pipe. The 
transverse momentum, of each particle is determined 
by its deflection angle in the azimuthal direction as mea- 
sured by the DC. The total particle momentum is recon- 
structed by projecting tracks back to the collision vertex 
through the magnetic field. The track reconstruction effi- 
ciency is approximately 98% and independent of pt with 
negligible centrality dependence. The particle momenta 
are measured with a resolution Sp/p = 0.007 O.Ollp, 
with p in GeV/c. The absolute momentum scale is cali- 
brated to 0.7% from the reconstructed proton mass using 
the PHENIX time-of-flight system. At high transverse 
momentum, a substantial background of electrons is pro- 
duced by photon conversion in the material between the 
beam pipe and drift chambers. To subtract this back- 
ground from photon conversion we use the ring imaging 
Cherenkov detector (RICH). The RICH is filled with CO2 
gas at atmospheric pressure and has a charged particle 
threshold 7th = 35 for emission of Cherenkov photons. 

To characterize the global parameters of the collision 
and its centrality, the PHENIX experiment uses beam- 
beam counters (BBCs) covering the pseudorapidity range 
3.0< |7y| < 3.9, zero degree calorimeters (ZDCs) and for- 
ward calorimeters (FCALs) located at \ri\ > 6. The lo- 
cations of these global detectors are schematically shown 
in Fig. [TJ Each BBC is an array of 64 Cherenkov coun- 
ters around the beam pipe and is positioned at 1.44 m 
upstream or downstream of the nominal vertex loca- 
tion. The information from the BBC is used for the 
event timing, vertex position, and centrality determina- 
tion. The two ZDCs are hadronic calorimeters which 
measure spectator neutrons They are located 18 m 
from the interaction point. At the top RHIC energy of 
100 GcV/nucleon, neutrons evaporated from the specta- 
tor remnants of the collision are emitted within 1 mrad 
from the colliding beam direction. Charged fragments 
and the noninteracted primary beam are bent by de- 
flecting (DX) magnets to much larger angles. The ZDC 
measures the total neutron energy within a small cone 
and with this provides the number of spectator neutrons 
from the interacting nucleus. The ZDC on the north side 
measured 100 GeV neutrons from dcuteron fragmenta- 
tion with a resolution of cr = 28 GeV. 

The forward calorimeters were installed before the 
d-l-Au run. The FCAL is a hadron calorimeter and con- 
sists of lead scintillating fiber modules originally used 
in BNL Alternating Gradient Synchotron (AGS) Exper- 
iment E864 [1^ rearranged into two 9 by 10 arrays. The 
only difference from the E864 experiment is the readout 
electronics, which are identical to PHENIX central arm 
electromagnetic calorimeter electronics. The size of each 
module is 10 x 10 x 117cm, the average tower density is 
9.6 g/cm'^, and the total length corresponds to 60 nuclear 
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FIG. 1: FCAL, ZDC, and BBC positions relative to the ver- 
tex. The vertical scale in the figure is arbitrary. The insert 
at the top shows the position of primary beam, spectator 
neutron and proton spots at the FCAL and ZDC locations; 
deuteron beam, spectator neutrons and protons go into the 
plane, the gold beam is coming out of the plane. 



interaction lengths. The two arrays are located 18 m from 
the interaction point along the beam pipes downstream 
of the first beam- line deflecting (DX) magnet. The DX 
magnets work as sweeping magnets for the spectator pro- 
tons. The FCAL measures the energy of the spectator 
protons. The resolution of the FCAL on the north side 
for measuring 100 GeV protons from the deuteron frag- 
mentation was cr = 40 GeV. 



B. Centrality analysis 

The present analysis is based on minimum bias events, 
defined by a coincidence of at least one photomultiplier 
each in the north and south BBCs. The data were taken 
for events with vertex position within |z| < 30cm along 
the beam axis. A total of 6.2 x 10^ events were analyzed, 
which corresponds to 1.6 nb^^ of total integrated lumi- 
nosity. At this vertex cut the minimum bias trigger cross 
section measured by BBC is 1.99 b ± 5.2% JJ,]. Thus 
at the trigger efficiency [12| of 88.5% ±4% we get the 
total inelastic d-|-Au cross section of 2.26 ± 0.1 b. Cen- 
trality classes in d+Au events were defined by charged 
particle multiplicity in the BBC south (BBCS) (the gold 
fragmentation side). We assume that the number of 
charged particles firing the BBCS is linearly proportional, 
on average, to the number of participants from the gold 
nucleus in the reaction. To check this, the BBCS re- 
sponse was simulated as a superposition of independent 



A'targ nucleon-nucleon type reactions, where iVtarg is the 
number of participating nuclcons in the struck gold nu- 
cleus. As a baseline for the BBCS modeling we use unbi- 
ased data from previous RHIC p + p runs where proton- 
proton inelastic collisions were selected with a trigger 
synchronized to fire whenever filled bunches crossed in 
the PHENIX interaction region. A'targ was calculated 
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using a Glauber model with a particular parame- 
ter set |21l] . The resulting distribution for minimum bias 




FIG. 2: (Color online) Comparison between the experimen- 
tal data hit distribution in BBCS in d-l-Au collisions (open 
circles) and the calculated BBCS hit distribution (solid line) . 



FIG. 4: Probability density distribution for the proton- 
neutron distance in the deuteron given by the square of the 
Hulthen wave function [2l| . 
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FIG. 3: Multiplicity distribution in BBCS, located on the gold 
nucleus fragmentation side. Four centrality classes for d-|-Au 
collisions are defined by slicing the BBCS distribution, shown 
with vertical lines. The same multiplicity cuts were used for 
the tagged sample of p-f Au and 7i-|-Au events, the summed 
distribution of which is shown in the lower histogram. 

The calculation appears to follow the general features 
of the data distribution, which supports our assumption 
that the number of BBCS firings is proportional to Nta^g- 
We define four centrality classes for d+Au collisions by 
slicing the BBCS multiplicity distribution into four re- 
gions, A, B, C, and D (Fig. [3]). These regions were se- 
lected to define in d-l-Au collisions four centrality classes 
0-20%, 20-40%, 40-60%, and 60-88%, respectively. We 
have taken into account the limited trigger efficiency for 



the peripheral collisions. The tagged sample, in which 
only one nucleon from the deuteron projectile interacts, 
has a differently shaped BBCS multiplicity distribution 
than the general data. Nevertheless, for both the tagged 
sample and the general data the same BBCS selections 
were used to define the collision centrality. 

In the Glauber model Monte Carlo simulation of d+Au 
collisions the deuteron wave function was represented by 
the Hulthen form ^^-^ The square of this wave function 
determines the probability distribution for the proton- 
neutron distance in the deuteron, as shown in Fig. 31 The 
deuteron is a large system with a mean proton-neutron 
distance of about 3 fm and a significant probability to 
be larger, fn the left side of Fig. O we present the cal- 
culated impact parameter distributions of rf+Au colli- 
sions and the tagged A^4-Au sample. On the right side 
of Fig. O the corresponding distributions of the number 
of collisions per participant nucleon from the interacting 
deuteron, v, are plotted. The parameter v is comparable 
to the number of collisions suffered by the proton in p+A 
experiments. The impact parameter b is defined as the 
distance between the centers of the colliding nuclei, Au 
and d. This means that for the tag sample, the inter- 
action point of participant nucleon from the deuteron is 



^ For the gold nucleus, we use the Woods-Saxon density distribu- 
tion with radius R=6.38 fm, diffuseness parameter a=0.54 fm 



and N-N cross section u]^^ 



42 mb. 

by a Hulthen wave function with a = 
fm-l, see Ref. [H. 



The deuteron is described 
0.228 fm-l and /3 = 1.18 
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FIG. 5: (Color online) Left: impact parameter distribution for the minimum bias d+Au collisions, for the most central events 
(centrality bin A) and for the tagged sample. For tagged events, the impact parameter was defined from the center of the 
deuteron. Right: the corresponding distribution of the number of collisions per participant nucleon from deuteron, v. 
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FIG. 6: (Color online) Scatter plot of ZDCN (vertical axis) 
and FCALN (horizontal axis) signals on the deuteron frag- 
mentation side. Solid lines show cuts which define the p+Au 
and n+Au collisions. 



closer to the Au center than the distance 6, and the spec- 
tator nucleon leaves the interaction region at a distance 
larger than h. We also did a Monte Carlo simulation with 
the most recent deuteron wave function parametrization 
from Ref. [23|. We found that the difference between 
the calculation with the Hulthen form and the most re- 
cent parametrization for the number of nucleon collisions 
does not exceed 1%. 

In Fig. [S] we illustrate our tag selection cuts. A p-t-Au 
collision event is tagged by the detection of a specta- 
tor neutron in the ZDC north (ZDCN) on the deuteron 
fragmentation side. Similarly, we use the FCAL north 
(FCALN) on the deuteron fragmentation side to detect a 
spectator proton and thereby tag n-|-Au collision events. 
The scatter plot in Fig. [6] shows the ZDCN and FCALN 
signals and has three distinct regions. Region 1 is defined 
as small or no signal in both the ZDCN and FCALN, 



which corresponds to the case in which both nucleons 
from the deuteron interact with the Au nucleus. Re- 
gion 2 has a small signal in the ZDCN and about 100 
GeV amplitude signal in the FCALN. This corresponds 
to tagged n-l-Au coUisions. Region 3 has a small signal 
in the FCALN and about 100 GeV energy release in the 
ZDCN. Events in this region are tagged p+A\x events. 
In Region 3, there is a small (anti) correlation between 
the ZDCN and FCALN. The reason for this is the close 
proximity of the FCALN to the ZDCN, see Fig. [H The 
ZDCN effectively acts as a secondary target for 100 GeV 
spectator neutrons. There is some contamination of sec- 
ondary particles produced in the ZDCN into the large 
volume FCALN. 

The purity of the tag samples from possible back- 
ground contamination was thoroughly investigated. Fig- 
ure [7] shows ZDCN and FCALN responses for the most 
central events. Background may contribute to the tag 
sample from the left side of the spectrum with low am- 
plitudes. To estimate the background contamination we 
fitted the left part of the spectrum using an exponen- 
tially falling function and fitted a Gaussian function to 
the detector response to the spectator nucleon. The thin 
smooth line in Fig. [7] shows the sum of the two func- 
tions. The background under the Gaussian peak was es- 
timated as the integral of the exponential function above 
our cut, which is marked by the vertical line and arrow. 
In central p-|-Au events we estimate 2.8% contamination 
in the ZDCN spectrum. For more peripheral events this 
contamination decreases and reaches 1.4% in the most 
peripheral bin. 

The background contamination in the FCALN defined 
n-|-Au sample is more complicated. For the most central 
events, as shown in Fig. [7] on the right, the background 
contributes 0.4% to the area above our cut. As centrality 
changes so too does the shape of the background spec- 
trum. Figure [H] illustrates the problem of background es- 
timation. Attempts to fit the background spectrum with 
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FIG. 7: ZDCN and FCALN spectra for the most central 
tagged events. Detector response to the spectator nucleon 
is fit with a Gauss function. Background from the left side 
of the spectra was fit to an exponential function and extrap- 
olated to the region above our cuts. All events to the right of 
the cuts, shown by vertical solid lines, were defined 
sample. 



an exponential function, as was done for central events, 
failed. The sum of the estimated background and the 
Gauss fit is well above the experimental data implying 
that the background spectrum falls much faster than a 
simple exponential function. However, the relevant ques- 
tion in assessing the background contamination is not its 
absolute or relative value, but rather how much it could 
distort our experimental data. The goal of this analysis is 
to measure the charged hadron spectrum in the PHENIX 
central arms. Thus, we can estimate the stability of our 
results measured by the central arms versus different cuts 
applied to the FCALN. We define the stability of the 
measurement, S, in a particular FCALN bin (Fig. [5]) as 
the ratio of the number of charged particle tracks per 
event with px > 0.5 GeV/c for this cut to a reference cut 
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FIG. 8: FCALN spectrum for the most peripheral tagged 
events, with the same line notation as Fig. [T] The plot illus- 
trates the difficulty of evaluating the background contribu- 
tion in the most peripheral collisions: background fit with an 
exponentially falling spectrum significantly overestimates the 
contribution to the events above our cut. 
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FIG. 9: (Color online) Stability of the results for different 
FCALN bins, see text for details. 



in which the signal to background ratio is the largest. 
For the reference cut we use 100 ± 20 GeV FCALN re- 
sponse (within region 3 as defined above). This is the 
most probable detector response for a spectator proton 
as seen in Fig. [S] We observe that in the region just to 
the left of our reference cut, where we may expect about 
equal amount of signal and background events, the pa- 
rameter S is only 3% smaller than unity, as seen in Fig.[Hl 
Within errors, S does not depend on track momentum. 
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Thus, in the worst case the background contamination ef- 
fect is less than 2% for events above our cut in the n+Au 
tagged sample. 

We do not correct the tagged sample results for any 
contamination. We average the two tagged samples to 
form a single nuclcon + Au (iV+Au) sample. 

For each of the four BBCS multiplicity regions shown 
in Fig. O using a Monte Carlo model within the Glauber 
multiple collision formalism [2^, we calculate for the 
d+Au and the tagged iV+Au collisions the average num- 
ber of participating nucleons iVpart, the nuclear overlap 
function Tab, the number of collisions A^coib and number 
of collisions per participant nucleon from the deuteron v. 
The calculated values are presented in Table [J The nu- 
clear overlap function Tab is defined as 

TAB{b) = J (fsTA{^TB{\b-s\), (1) 

where the integration is performed over the element of 
overlapping area cP s, and s — {x, y) is a, vector in the 
transverse plane of interacting nuclei at the impact pa- 
rameter b between the centers of the nuclei. For nucleus 
A the nuclear thickness function TA{b) is defined as 

TAib) - [ dzpA{b,z). (2) 



TABLE I: Total number of participants A'part , number of colli- 
sions A'coii, nuclear overlap function Tab, see Eq. 1101 average 
number of collisions per participant nucleon from deuteron 
V, and the BBC bin correction factor for different centrality 
classes. 



Cent, bin 




(A'coii) 


{Tab), mb-i 


u 


Obbc 


A 


15. Oil. 


15.4±1.0 


0.37±0.02 


7.5±0.5 


0.95±0.03 


B 


10.4±0.4 


10.6±0.7 


0.25±0.02 


5.6±0.4 


0.99±0.01 


C 


7.0±0.6 


7.0±0.6 


0.17±0.01 


4.0±0.3 


1.03±0.01 


D 


3.2±0.3 


3.1±0.3 


0.07±0.01 


2.2±0.2 


1.04±0.03 


tag A 


10.6±0.7 


9.6±0.7 


0.23±0.02 


9.6±0.7 


0.93±0.03 


tag B 


8.0±0.6 


7.0±0.6 


0.17±0.02 


7.0±0.6 


0.95±0.02 


tag C 


5.6±0.3 


4.6 ±0.3 


O.llifcO.Ol 


4.6±0.3 


0.95±0.02 


tag D 


3.1±0.2 


2.1±0.2 


0.05±0.01 


2.1±0.2 


0.97±0.04 



Normalization of Tab (b) is done by integration over all 
impact parameters: 

j SbTAB{b)^AB. (3) 

The average number of binary inelastic nucleon-nucleon 
collisions at impact parameter b was calculated from 
TAB{b) as 

{Ncoll) =<yNNTAB{b), (4) 

where a^M is the inelastic nucleon-nucleon cross sec- 
tion [m. 

An additional multiplicative correction factor, CbbCi 
see TablelH has been applied to the data for different cen- 
trality bin selections [23|. This correction addresses two 



effects, each of which distorts the centrality classification 
in the opposite direction. Because of natural fluctuations 
in the number of produced charged particles at a partic- 
ular impact parameter, the BBCS centrality selections 
have imperfect resolution. In the case of a steeply falling 
BBCS multiplicity spectrum, especially for the tagged 
sample (see Fig. [3]), there is a contamination of periph- 
eral collision events into a more central event class. This 
effectively decreases the actual number of A'^coii and par- 
ticle production in the central events. The second effect 
is due to the BBC coincidence required for a PHENIX 
event trigger. In all calculations we used for the inelastic 
nucleon-nucleon cross section gnn = 42 mb. Actually, 
this cross section has three components: nondiffractive 
collisions with 28 mb, single diffractive collisions with 
10 mb, and double diffractive collisions with 4 mb cross 
section [24]. From Monte Carlo simulation of p+p col- 
lisions we found BBC trigger efficiencies of about (72 
± 1)%, (7 ± 1)% and (32 ± 1)%, for each process re- 
spectively. In the Monte Carlo simulation we use the 
PYTHIA 6.2 event generator [13] and the GEANT simu- 
lation of the BBC detector. Single and double diffractive 
collisions produce particles dominantly near the beam 
rapidity and thus have a small probability for particle 
production in the BBC acceptance of 3.0< \ri\ < 3.9, and 
an even smaller probability at midrapidity. Therefore, 
the BBC trigger is biased to the nondiffractive collisions, 
which have larger particle production at midrapidity. In 
more central events as the number of nucleon-nucleon 
collisions increases the probability that there will be at 
least one nondiffractive collision approaches 100%; thus 
the BBC bias becomes negligible in central events. 



C. Charged hadron analysis 

The present analysis follows the methods of the anal- 
ysis described in Ref. The majority of background 
tracks are particles with low momenta, which in traveling 
from DC to PC3 undergo multiple scattering and are ad- 
ditionally deflected in the residual magnetic field behind 
the DC. 

To minimize this background we employ a track match- 
ing cut in the PC3 that rejects tracks whose displacement 
in the </) or z direction, and , respectively, is greater 
than 2.5 standard deviations. In addition we make a fidu- 
cial cut around the z vertex determined by the BBC. De- 
spite these veto cuts, there is still significant background 
contamination for above 4 GeV/cwhich must be sub- 
tracted. The main sources of the remaining background 
are e~^e^ pairs from photon conversions in the material 
between the collision vertex and the drift chamber, and 
secondary particles produced by hadron decays. 

To distinguish these backgrounds, we use the RICH de- 
tector to divide all tracks into two subsets: tracks with an 
associated RICH signal, Njf , and tracks with no signal in 
the RICH, Npfji. Tracks with at least one hit in the RICH 
contain high pious and conversion electrons. For re- 
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FIG. 10: (Color online) PCS Dtf, distributions for the conversion subtraction (left) and the decay subtraction (right) shown for 
minimum bias with 6 < pr < 7 GeV/c. The raw distributions, Nr and Nnr for conversion and decay respectively, are shown 
with dotted lines. The estimated conversion background shown with a dashed line is N^/ Re- The estimated decay background, 
shown with a dashed line, was obtained by scaling the PCS distribution of Nnr tracks with pr > 10.5 GeV/cbased on the 
(4(T < \D^\ < 10a) region. In both plots, the signal is the raw minus the estimated background distributions. The vertical bars 
show the track matching cuts at ± 2.5 D^. 



constructed electrons with momentum above 150MeV/c, 
the average number of pliotomultiplier tube (PMT) hits 
in the RICH associated vifith the track is (iVpMT) ~ 4.5. 

The RICH detects more than 99% of all conversion 
electrons for A'pmt > 1- At this threshold the RICH 
also detects pions with pt> 4.8GeV/c, but the number 
of associated PMTs for pions reaches its asymptotic value 
wen above 10 GeV/c; for a 10 GeV/c pion (A^pmt) = 3.6. 
Therefore we label tracks with TVpmt > 5 as electron 
tracks, N^, which compose some fraction, R^, of conver- 
sion electrons. To calculate this fraction, we take ad- 
vantage of the deflection of conversion electrons in the 
magnetic field between the DC and PCS. This deflection 
leads to poor track matching in PCS which distinguishes 
electrons from true high pions. We define poor PCS 
track matching as a displacement of more than four stan- 
dard deviations in the (j) direction. Thus we measure the 
value of Re as the fraction of Nfi with 4ct < jlJ^I < lOcr. 
For minimum bias events we find Re = 0.41 ± 0.01. The 
real pion signal Sr in the Nn sample, is calculated for 
each pt bin as 

Sr = Nh-^ (5) 

He 

The PCS distribution for the conversion subtraction is 
shown on the left side of Fig. (TUl The conversion subtrac- 
tion is performed independently in each centrality bin. 
The definition of Ng (tracks with A'pmt > 5) does not 
perfectly select electron tracks, as some fraction of pions 
satisfies the cut. This leads to a fraction of authentic pi- 
ons, which have A^pmt > 5, being subtracted along with 
the conversion electrons. This fraction is small below 
7GeV/c, but increases rapidly for higher px- We calcu- 
late a correction factor to address this over-subtraction, 



using a Monte Carlo simulation of the detector. The un- 
certainty associated with this correction, loss , is shown 
later in Table [ml 
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FIG. 11: Feed-down contamination as the ratio of protons 
and charged pions from feed-down decays to total detected 
charged hadrons. The low px points shown with triangles are 
calculated using the fraction of p + K particles to hadrons 
measured in Ref. [2^. The higher pT points are calculated 
using the various ratios of Eq. [S] The solid line is a fit to both 
sets of points and is used as Cfcod. It is bracketed by dashed 
lines showing the assigned 50% systematic uncertainty. 

The A^-I-Au sample is averaged over p-fAu and n+Au 
events. This is based on the good agreement between the 
yields of the two tagged samples, as shown in Fig. [T31 

The hadron decay background is of two types: "decay" 
and "feed down." The decay background is produced by 
TT and K decays far from the source and thus with re- 
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FIG. 12: (Color online) Charged hadron spectra for centrality selected d+Au collisions and A'^+Au collisions. The error bars 
represent statistical uncertainties only, although they are generally smaller than the data point symbols. Points from different 
centrality bins are scaled sequentially by a factor of 10. (See Table |T] for centrality class definitions.) 



constructed momenta different from their true momenta, 
whereas the feed-down background is produced by weak 
decays of short-lived particles, mostly and A parti- 
cles near the event vertex with apparent momenta close 
to their true momenta. 

From tracks with no RICH signal, the N^n sam- 
ple, we define a subsample by selecting tracks with px 
> 10.5GeV/c, a px region which is almost exclusively 
background. We expect that its shape in will be the 
same as the background in the lower px region. Within 
this subsample we calculate the ratio, i?docay, of tracks 
which pass the PCS cut {\D^\ < 2.5a) to those with a 
poor match (4(7 < \D^\ < 10a). For minimum bias events 
Ji^docay = 0.55 ± 0.03. For each momentum bin, the total 
decay background is then obtained by multiplying the 
Nnr tracks with poor PC3 matching {Aa < \D^\ < 10a) 
by i?dccay The PCS distribution for the decay subtrac- 
tion is shown on the right side of Fig. [TOl The decay 
background as a function of px is measured and sub- 
tracted independently in each centrality bin. 

The feed-down subtraction addresses the detected tt 
and p particles that were produced in the decays of 



and A particles, so we define the total feed-down contam- 
ination as: 



(p + tt) 



feed 



^detected 



(6) 



averaging over charged particle and antiparticle yields. 
For feed-down estimation, we have no statistical recourse 
and must resort to simulation to find the contamination. 
We assume that the spectral shapes of the K^J and A fol- 
low the shapes of the charged kaon and proton spectra, 
respectively. There is good agreement of Monte Carlo 
simulations in d-l-Au [2a | and minimum bias Au-|-Au Q 
for proton feed down. We therefore use the Au-|-Au 
Monte Carlo simulation, which also includes to pion 
processes, to obtain the ratio (p-|-7r)^°°'^/(p-|-iir )'^®*®'^*®'^ = 
0.2, after the decay background subtraction. To make use 
of this ratio we rewrite the contamination as: 



detected 



Cfeed — 



{p + tt)"'''^ (p + K) 

{j) -|- J^^dctcctcd ^detected 



(7) 



To find the contamination we use the fraction oi p + K 
particles in our measured hadron sample. To do so, for 
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Pt less than 2.5GeV/c we use the PHENIX published 
data on identified hadron production in d+Au ^26J. We 
exphcitly calculate the (p+ii:)dotoctcd//jdotoctcd ^.j^^-Jq ^^^^ 

those data. At higher assuming tt" has the same yield 
as TT^, we subdivide this ratio as 



K 



h 



Kit 



K 



~h 



(8) 



so that the right-hand side consists of all measured quan- 
tities. The p/tt ratio is taken from STAR measure- 
ments [131 scaled to match the PHENIX data [26| in their 
common px re gion . The K/i:^ ratio is calculated from 
Refs. [l^ and [2^. The ratio of tt^s to hadrons is calcu- 
lated from the charged hadron measurement of this anal- 
ysis and 7r° measurements in Ref. [12] • From both the 
low and high regions Cfeod is calculated to be 9.6%. 
To this factor we assign a 50% systematic uncertainty 
based on uncertainty of the various particle ratios and 
the Monte Carlo simulation. Cfccd is shown in Fig. [TT] 
We correct for the effects of feed-down decay by multiply- 
ing the spectra remaining after the conversion and decay 
subtraction by 1 — Cfccd- 

Following the background subtraction, we constructed 
a single, px dependent correction function to correct 
the hadron spectra for acceptance, decay in flight, re- 
construction efficiency, and momentum resolution. The 
correction is determined by using a Monte Carlo sim- 
ulation pJi] of the PHENIX detector. The correction 
function is necessarily particle species dependent to take 
into account multiple scattering and decays, therefore 
we calculate separate correction functions for tt"*", tt", 
K~ , and p~ . The individual functions are 
weighted by the particle px spectra measured in periph- 
eral Au-|-Au collisions [2^| to form a single correction 
factor, Cmc(pt), subject to a systematic uncertainty, 
(^MCweightjStemming from uncertainty on the particle 
mixture [l4|- For absolute normalization of the spec- 
tra we match the geometrical acceptance of the Monte 
Carlo simulation with the actual acceptance of the data. 
To obtain the charged hadron yield we multiply the 
background-subtracted spectra by the correction func- 
tion. We normalize each centrality bin by dividing by 
the number of events, and each momentum bin by di- 
viding by its bin width. Each data point is corrected so 
its value corresponds to the bin center. We define the 
invariant hadron yield as 



1 



Ncvt 27rpT dpT d-q 



d^N 



bkg— subtracted 



2-KpT dpT drj ^ 
(1 — Cfocd) Cmc^t) Cbbc 



(9) 



All of the preceding steps are applied to the tagged 
N+Au (nucleon -I- gold) sample as well as to the gen- 
eral d-l-Au sample. 

We further examine the d+Au invariant yields by com- 
paring them with the invariant yield from p+p collisions. 
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FIG. 13: (Color online) Ratio of n-|~Au to p-|-Au invariant 
yield per event shown in the most central and peripheral bins. 
The error bars are statistical only, as all systematic uncer- 
tainties cancel in the ratio. (See Table |I] for centrality class 
definitions.) 



Previous experiments have demonstrated suppression of 
charged hadron and pion yields in Au-|-Au collisions [^; 
such suppression can be quantified by the nuclear mod- 
ification factor Rab- For any collision of nuclei A + B, 
Rab is calculated for each centrality class as the ratio of 
the yield in A + B collision to the cross section in p + p 
collisions scaled by the nuclear overlap function (Tab), 
that is. 



Rab{pt) = 



(l/TVevt) d^N^+^/dpT dy 
{Tab) d^aP+P/dpT dr; 



(10) 



{Tab) is determined by the density distribution in the 
nuclei A and B and is averaged over the impact parame- 
ter range within a particular centrality class. In our case, 
nucleus A refers to the deuteron (or the single nucleon 
in tagged events), and nucleus B refers to the gold nu- 
cleus. {Tab) values are presented in Table |T1 Using the 
PHENIX p + p spectra [29] , we thus calculate -R^Au and 

RnAvl- 



TABLE II: Systematic uncertainties that are constant for all 
Pt- 



Source 



Unccrtainty(%) 



Geometric acceptance correction 
Track matching 
Run-by-run variation 
Feed-down correction 

^MCwcight 

Total 



2.9 
2.2 
5 

4.8 
3.7 
8.7 



We estimate systematic uncertainties in the methods 
and assumptions of our analysis as displayed in Tables |TT] 
and IIIII Table [Til shows uncertainties in the spectra that 
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FIG. 14: (Color online) From left to right: RdAu, RnAu as functions of pr- The bars represent statistical uncertainty (often 
smaller than the data point symbols), the shaded boxes on each point represent systematic uncertainties that change with pr, 
and the black boxes on the left show systematic uncertainties that do not change with pr- (See Table |I] for centrality class 
definitions.) 



TABLE III: Systematic uncertainties that vary with pr- 
Background subtraction (Bckg. sub.) uncertainties refers to 
the 0-88% d+Au spectra; the uncertainties are greater in the 
more peripheral d+Au and more central 7V+Au spectra. 



PT (GeV/c) Mom. Mom. Bckg. S^i^ss (%) Total(%) 
res. (%) scale (%) sub. (%) 



<4.5 


<0.5 


<3.2 


<0.1 


<0.3 


<3.3 


4.5-5.5 


<0.6 


3.3 


0.5 


0.5 


3.4 


5.5-6.5 


0.8 


3.5 


1.4 


1.1 


4.0 


6.5-7.5 


1.0 


3.6 


2.0 


3.6 


5.6 


7.5-8.5 


1.4 


3.7 


4.9 


6.9 


9.3 


8.5-9.5 


1.8 


3.8 


11.9 


13.9 


18.8 



although the uncertainty from the feed-down correction 
partially cancels in the Rab because of an analogous un- 
certainty in the p + p reference spectrum. The back- 
ground subtraction and momentum scale uncertainties 
only partially cancel in the Rab] the momentum res- 
olution uncertainty and (5^ioss fully cancel. The p + p 
reference spectrum also introduces an uncertainty that 
ranges from 10.7% in the low px bins to 11.3% at > 2 
GeV/c. The 7r° analysis [13| we used in the hadron/7r° 
ratios applies the same BBC bias correction factors as the 
present analysis and therefore the uncertainty stemming 
from this correction factor is canceled. 



III. RESULTS 



do not vary with px, and Table HIT] shows uncertainties 
that vary with px- All the uncertainties listed in Ta- 
ble [IT] propagate into the Rab and the hadron/7r" ratio, 



The fully corrected px distributions of (ft,+ + h^)/2 for 
d-|-Au and 7V-|-Au collisions for minimum bias and four 
centrality classes are shown in Fig. [T^l The d-f Au data 
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FIG. 15: Ccolor online) -R^Au (open symbols) and RnAu 
(closed symbols) values averaged in three momentum ranges 
as functions of v—l. Point-by-point uncertainty bars show the 
quadratic sum of the statistical uncertainty in Rab and the 
systematic uncertainties that change with v. The box on the 
far left represents the size of the systematic uncertainty that 
does not change with v. Horizontal bars show the uncertainty 
in the value of v for each centrality class. 



are in good agreement with the sum of identified hadrons 
published in Ref. [2G]. 

Figure [HI shows RdAu and Rn Au- As expected, in the 
most peripheral iV-|-Au bin, with A^coii = 2.1, Rn Au is 
consistent with unity. i?c;Au and Rn Au are in agreement 
within our uncertainty bounds. 

The enhancement of the hadron yield relative to p + p 
collisions has previously been observed in lower energy 
p+A collisions [l3, [HI, and is called the Cronin ef- 
fect. We observe that RdAu and Rn Au are system- 
atically larger than unity in the momentum range be- 
tween 1.5 GeV/c and 5 GeV/c with maximum amplitude 
around 1.3. 

There are many theoretical models with very different 
assumptions about initial state effects, which describe the 
Cronin effect [sol . [sil . [33 , [ssj . All the models agree that 
there is at least one additional scattering of the initial 
nucleon or parton while propagating though the target 
nucleus. This scattering increases the intrinsic transverse 
momentum of the colliding parton, and leads to a broad- 
ening of the parton distribution. We can parametrize 
the effect of this broadening by writing the mean value 
of parton intrinsic momentum kx as 



{kl)A: 



(11) 



where {k'^)pp is the square of the initial parton trans- 
verse momentum in the proton, {k'^)A is an additional 
momentum squared after rescattering, and {k^)pA is the 
final broadened width. Most of the models differ on the 
assumption they use to describe (fc|,)^: whether there is 
a single hard scattering (3^ or a sum of small sequential 
rescatterings |3l| that produces the additional kx- 
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FIG. 16: (Color online) (/i^ + /i~)/27r'' ratios as functions of 
transverse momentum from d-|-Au collisions in four centrality 
bins. Open circles are the charged hadron spectra from the 
present analysis divided by 7r° data from Ref. [lH. Bars indi- 
cate statistical uncertainties, and the shaded boxes systematic 
uncertainties. Triangles are the (p* -I- + 7r*)/7r* ratios 
from Ref. [2^, with statistical uncertainties roughly the size 
of the symbols. (See Table |I] for centrality class definitions.) 



Common to the models is that {k'^)A is a function of 
the number of sequential nucleon-nucleon collisions, v. 
For impact parameter 6, {k"^) a can be written as: 



(4)^(6) = H{v{b) - 1) 



(12) 



where H is the square of the average momentum acquired 
'mv ~ 1 rescatterings. For a single hard scattering model, 
{k'^)A should saturate at = 2. We therefore investigate 
the shape of RdAu as a function of v to illuminate the 
underlying process. The centrality selection of our data 
and the tagged 7V-|-Au sample allow us to investigate 
precisely the effect of the collision geometry. We use 
1^ = (A^coii/A^^a™*"") (in ^^+Au collisions = (N.^n)) 
to look explicitly at the impact parameter dependence of 
the nuclear modification factor. 

In Fig. I15[ we plot -R^Au and Rn Au as a function of 
(i/ — 1). The values of ly are presented in Table HI 

Three transverse momentum regions were selected to 
study the dependence of Rab on v. 2.8< pt < 6.0, 1.5< 
Pt < 2.7, and 0.6< px < 1-0 GeV/c. In the low pT re- 
gion, we expect scaling with the number of participating 
nucleons rather than with the number of binary colli- 
sions; therefore Rab is less than unity. The Cronin ef- 
fect is observed in the 2.8< pt < 6.0 GeV/c region, where 



14 



within the hmits of our uncertainties, it is independent 
of the number of additional scatterings {ly ~ 1). In the 
intermediate pt region, the data show httle to no Cronin 
enhancement, thus confirming scahng with the number 
of binary colhsions. Just as RdAuipT) and RnAuIpt) 
matched very closely, so do RdAui'^) and RNAui'^)- 

As discussed above, previous experiments found a 
larger Cronin enhancement for protons than pions [ll[. 
Recently the PHENIX Collaboration published identified 
pion, kaon, and proton production data from ^^snn ~ 
200 GeV d+Au collisions [1^ , in which a similar Cronin 
behavior was observed. RdAu for protons reaches about 
1.8 at Pt ~ 3GeV/c, whereas pion RdAu is measured to 
be about 1.1 for transverse momentum between 2 and 
2.6GeV/c. 

At the time these data were taken, the PHENIX exper- 
iment did not have the capability to measure identified 
charged pions and protons at momenta above 2.6GeV/c 
and 4 GeV/c, respectively. To extend the trends observed 
in these low momentum regions [1^ , we calculate the ra- 
tio of the charged hadrons measured in the present analy- 
sis to the 7r° spectra from Ref. [l^ , alongside the ratio of 
{p^ + +tt^)/tt^ from [2^. These ratios are presented 
in Fig. [m 

In the transverse momentum region common to the 
two presented ratios there is strong agreement between 
the analyses. The -I- h~)/2Tr'~' ratio is independent 
of Pt above 2GeV/c where the identified particle data 
end, implying that the RdAu and particle ratio trends ob- 
served at low transverse momentum continue at higher 
transverse momentum. The average value of 1.58 ± 0.03 
of the -I- h~)/2TT° ratio for px above 2GeV/c in the 
peripheral D centrality bin agrees well with the value of 
1.59 obtained from lower energy collisions in the GERN 
Intersecting Storage Rings (ISR) ^3^]. As found in Ref. 
Q for Au-|-Au collisions, this value rises for more central 
events: we find average values of 1.78 ± 0.02, 1.77± 0.03, 
and 1.72 ± 0.03 for centrality bins A, B, and G, respec- 
tively. There is an additional 11% systematic uncertainty 
common to all four values. The centrality dependence 
implies moderate medium modification effects in central 
^/SNN — 200 GeV collisions, even in the d+Au system, 
that increase the production of protons and kaons rel- 
ative to pions. The PHENK measurement of particle 
species dependent -R^Au in [26| at lower p^ suggests that 
the increased particle production relative to pions is dom- 
inantly proton and not kaon production. 

IV. SUMMARY 

We have measured the centrality dependence of 
charged hadron production at midrapidity in d-l-Au col- 
lisions, as well as in p+Au and n-|-Au collisions that are 
tagged by a spectator nucleon from the deuteron nucleus. 
The hadron yields in p+Au and n+Au collisions are iden- 



tical within our experimental uncertainty. Using the p+p 
data from same energy, we calculated the nuclear mod- 
ification factors RdAu and RnAu for various centrality 
selections. Within this analysis and its experimental un- 
certainty, there is no difference between RdAu and RnAu 
. Instead of a strong suppression as predicted by some 
color glass condensate models an excess of hadron 
production is seen at >2 GeV/c, consistent with en- 
hancement due to the Cronin effect. The magnitude of 
the Cronin effect is independent of the number of ad- 
ditional scatterings {v - 1) within the limits of our un- 
certainties. We also studied the ratio of charged hadron 
yield to pion yield. Wc found that the charged pions ac- 
count for about 60% of the charged hadrons at >3 
GeV/c, with a slightly larger value in central d+Au col- 
lisions. This implies that RdAu for protons and kaons 
remains close to one at higher px- 
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